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Abstract 



This work is related to different questions within cosmology. The principal idea 
herein is to develop cosmological knowledge making use of the analyses of obser- 
vational data in order to find the values of the matter density and vacuum 
energy density Q\. 

Data fitting is carried out using two statistical methods, maximum 
likelihood. The data analysis exhibits that a low density and flat Universe is 
strongly favoured. 

Applying the value found for clusters of galaxies, we demonstrate that 
clusters have very little room for baryonic dark matter. An upper limit to the 
small but non-negligible sum of baryonic dark matter and galaxy mass can be 
estimated, requiring the use of special statistics. 

A Toroidal Black Hole (TBH) study, in contrast to the Spherical Black Hole 
(SBH), shows that the TBH can be used as an important tool in explaining AGN 
phenomena. 
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Chapter 1 

Cosmological Models and 
Parameters 



Measurement of global cosmological parameters is one of the key challenges in 
current cosmology research. People are interested in knowing in what kind of 
Universe we live. Only a precise measurement of density parameters can tell us 
the critical nature of the Universe. The parameters determined from the different 
observations are strongly model dependent and also contradictory to each other. 
Resolution of contradiction is another big challenge to the scientist. A lot of 
interesting discussions is reported in this thesis, based on recent output from the 
cosmological research. 

According to the cosmological principle, the Universe is homogeneous and 
isotropic. Most of the cosmological models are based on Friedman equations, 
describing a Universe in expansion or contraction. Einstein sought a solution in 
which the Universe would be static and eternal because the Universe known at 
that time comprised only the Milky Way, which clearly was not expanding or 
contracting. In this Chapter a few models will be discussed that are remarkable 
in modern cosmology. 

This thesis deals with Cosmological Parameters and Black Holes. Gravita- 
tional matter is composed of baryonic matter, baryonic dark matter and dark 
matter of unknown composition. The matter density (0,^), baryonic density (f^b), 
Hubble constant (Hq) and the age of the Universe (to) will also be discussed. A 
discussion on the cosmological constant (A) and its vacuum energy density (Oa) 
will be presented in the next Chapter. A simple picture of different cosmologies 



are portrayed in Fig. 1.1. Some aspects of Black Hole physics will be discussed 



1 



in Chapter 6. 

1.1 Einstein's Model 

Modern cosmological models began with Einstein's static cosmology. Einstein 
proposed that the gravity is not a force but merely a manifestation of free mo- 
tion in space-time. In addition, the geometry of space-time is determined by the 
energy in the Universe. The energy is associated with mass, radiation and pres- 
sure. A cosmological constant A was introduced by Einstein into the equation 
of General Relativity to allow for a stationary solution, and it was introduced 
before the concept of Big Bang. However, the idea of a static Universe fails to 
agree with observations, because the Universe is expanding linearly with time. 

1.2 The de Sitter Model 

In 1917, Willem de Sitter discovered a cosmological model differing from Ein- 
stein's. He observed that the light from distant objects becomes redder as the 
distance increases. There is no matter in this model. The space-time is flat, 
and the space expands exponentially because of the A force. This is one of the 
earliest models that was considered when the modern science of cosmology was 
in its infancy. 

1.3 The Einstein-de Sitter Model 

According to Einstein and de Sitter, the cosmological constant should be set equal 
to zero, and they derived a homogeneous and isotropic model. They assumed that 
the spatial curvature of the Universe is neither positive nor negative but rather 
zero. The special geometry of the Einstein-de Sitter Universe is Euclidean, but 
the space-time is not globally flat. People with a philosophical bent have long 
considered it as the most fitting candidate to describe the actual Universe. Strong 
theoretical support for this viewpoint came from particle physics, but as yet not 
definitive, astronomical observations also supported this model. 

1.4 The Friedman- Lemaitre Model 

The standard model of cosmology describes the Universe as expanding at present. 
The formulation and prediction of a Big Bang expansion for the Universe is re- 
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markable. Both Friedman and Lemaitre in different years, independently dis- 
covered the solution to Einstein's equations of gravitation which described an 
expanding Universe. The expansion could either continue forever, or eventually 
reverse into a phase of contraction. 

However, the Fricdman-Lemaitre evolutionary model could not predict explic- 
itly the curvature of the expansion, permitting three possible solutions: either 
hyperbolic, flat or spherical expansion. Friedman's equations may be written 
(Kolb k Turner 1990, Roos 1997, Bergstrom & Goober 1999) as 

i?2 + kc^ A 87rG/9 
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(1.1) 



2R ^ B? + kc? ^ _ SnGp 

B? c2 ' ^-^-^^ 

where R is the scale factor of the Universe, k is the curvature parameter, A is the 
cosmological constant, G is the Newtonian constant, p is the total energy density 
and p is the total pressure. 

There is a relation between pressure and density, called equation of state 

p = a-p (1.3) 

where a is a constant, different in different eras of the Universe. In the radiation 
dominated era P ~ and in the matter dominated era P For the de 

Sitter model a = -1, but for interacting fields and topological defects a can vary 
between -1 and 0. Since recent supernovae observations show that the Universe 
is accelerating, one has -1 < a < -\. Models with a in this range some people 
call Quintessence models. 

The critical density of the Universe at redshift z is 

t ^ 3//2(^) 

= -MG- ^'-'^ 
where H is the Hubble constant. The vacuum energy density parameter Cl\ is 
related to the cosmological constant A by 

= A = ^, (1-5) 
6Hq Pc 

where the notation zero means the value at present. The matter density param- 
eter of the Universe is 

nm = ^. (1.6) 

Pc 
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The total density parameter of the Universe is the sum of matter density 
parameter and the vacuum energy density parameter 



nO = ^m + ^A- (1.7) 

The flat condition of the Universe is [/;; = in Eqs. (1.1) and ( |1.2| )] 

n^ + nA = i (1.8) 

We can arrive at a relation among ffoi^mj and tQ, from the Friedman- 
Lemaitre model, as 

1 

to = yr dx[{l-Qm-^A)+^mX-^ + ^AXY^^^. (1.9) 
-no Jo 

This equation will be used to determine the age of the Universe. There is no 
analytical solution to this integration but it can easily be evaluated numerically. 
To determine to, we need to have precise values of other parameters. The eval- 
uation of ^Im, ^A can be found in Papers (II), (III), (IV) and also in the next 
Chapters. 



1.5 Density Parameters 

Generally, matter density and vacuum energy density are treated as Universal 
density parameters, and the total matter density is the combination of different 
kinds of matter. So far as we know, there are two kinds of matter in the Universe, 
baryonic matter and dark matter. Thus the total mass density parameter r^rn is 
the sum of these two. The candidates of dark matter and their nature will be 
discussed in the Chapter 5. For a long time, the physicists had neglected Qa, 
simply it was considered that = and ^Im = 1. The situation has changed 
in the last few years by the enormous data collection of powerful telescopes and 
by the success of satellite missions. Observations suggest that this is not our real 
Universe. 

Due to the expansion of the Universe the mass density is decreasing and the 
vacuum energy density is now dominating. The vacuum energy is related to the 
cosmological constant, Eq. (|1.5| ). Most recent measurements agree that the value 
of present matter density $7^ ~ | in the flat case. It is the proper time for 
the astrophysicists and cosmologists to have a direct constraint on the density 
parameters. The constraints on matter density and vacuum energy density 
JIa are extensively described in the next Chapters. 
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1.6 Baryonic Density 



The value of the universal baryonic density parameter 0^ follows from standard 
Big Bang Nucleosynthesis (BBN) arguments, and from the observed abundances 
of ^He, D, ^He, and ''Li (cf. eg. Sarkar 1999), in particular from the low deuterium 
abundance measured by Buries et al. (1999, 2001). Their estimate is 

Qbh"^ = 0.020 ± 0.001 (68% CL) (1.10) 

where Hq = lOO/i. Thus, taking this central value and the Hubble constant from 
Gibson &: Brook (2001) we get today's universal baryonic density ~ 4%. The 
BBN value is in disagreement with the recent observations of CMB anisotropy 
which yields (Jaffe et al. 2001) 

Qbh^ = Om^tofol (68% CL). (1.11) 

Note that BBN provides a probe of the Universal abundance of baryons when 
the Universe was only a few minutes old. Observations of CMB anisotropy probe 
the baryon abundance when the Universe was 3-5 hundred thousand years old, 
and SNe la supernovae and clusters of galaxies observations probe a more recent 
past, when the Universe was several billion years old. 

Most recently, the Degree Angular Scale Interferometer (DASI) has measured 
the angular power spectrum of the Cosmic Microwave Background anisotropy 
over the range 100 < I > 900 (Pryke et al. 2001). Here, the second peak is 
more pronounced than found BOOMERANG (de Bernardis et al. 2000) and 
MAXIMA-1 (Balbi et al. 2000), and the contradiction with BBN can be resolved. 
DASI (Pryke et al. 2001) gives a new value of baryonic density as 

Qbh^ = 0.0221°:°°^ (68% CL), (1.12) 

which is a good agreement with BBN (Buries et al. 2001). DASI has also inde- 
pendently determined Q^n = 0.40±0.15 and 17a = 0.60±0.15 (68% CL), although 
the errors ber is large and no contour plot is given. Adding the data at higher 
I BOOMERANG (de Bernardis 2001) has reported the value of rj^/i^, their new 



value is exactly same as (DASI 2001) Eq. ( |1.12D . 

Assuming that the BBN value for Qbh^ in Eq. dTToD (Buries et al. 2001) is 
correct, and taking the value of the Hubble constant h = 0.73 ± 0.05 Gibson &; 
Brook (2001), one can obtains a constraint on the baryonic density parameters, 

= 0.038 ±0.002. 
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1.7 The Hubble Constant 



In 1929, Edwin Hubble discovered that the distant galaxies are moving away from 
us. A simple mathematical relation between the Hubble constant and distant 
objects is Hq = v/d, where v is the galaxy's radial velocity and d is the galaxy's 
distance from the earth. This is one of the important parameter because it is 
directly related to the age of the Universe, and it is also related to the other 
cosmological parameters. The Hubble constant is well determined from the HST 
Key Project on the Extragalactic Distance Scale (Mould et al. 2000), obtained 
by the combination of four methods (SBF, FP, Tully-Fisher, and SNe la). Their 
combined result is Hq = 71 ib 6 km s~^Mpc^^. A new value of the Hubble 
constant is presented by Gibson & Brook (2001) performing the re-analysis of the 
combined Calan-Tololo and Center for Astrophysics (CfA) type SNe la datasets. 
The fit is extremely good with 3% statistical and 10% systematic errors. With a 
calibration to the corrected Hubble Diagrams by seven high-quality nearby SNe 
la the Hubble constant becomes 



i?o = 73 ± 2(r) ± 7{s) km s"^Mpc"^ . 



(1.13) 



However, the results agree with each other at la level. The Hubble constant 



as measured by different methods summarized in table 1.1 



Method 







Error (%) References 



Type la super novae 73 

Combined HST methods 71 

Tully-Fisher clusters 71 

FP clusters 82 

SBF clusters 70 

Type II super novae 72 
Cepheids, Metallicity-corrected 68*) 

S-Z effect 66 



±2(r) ± 7(s) 
±2(r) ±6(s) 
±3(r) ± 7{s) 
±6(r) ±9(s) 
±5(r) ±6(s) 
±9(r) ± 7{s) 
±5(r&s) 
±13(r) ± 15(s) 



ii 



PI 



111] 



11 



m 



m 



Table 1.1: The Hubble constant in units of km s~^Mpc~^. Random (r) and 
systematic (s) errors are given at la confidence level. *) Taking into account 
more recent calibrations of the LMC distance, this value increases to 74 (cf. eg. 
Gibson & Brook 2001). 
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1.8 Age of the Universe 



The absolute age of the Universe to depends on {h,Qm,^A)- The mass den- 
sity of the Universe and the cosmological constant are known from the 
present analysis. Now we can apply a constraint on to, since our analysis strongly 
preferred a flat geometry in the Universe and this constraint is valid is the flat 
case also. In Papers II & III we have determined the age of the Universe to 
be to = 13.5 lb 1.3 (0.68/h) Gyr, where we use the earlier Hubble constant 
Ho = 68 ±5 km s^^Mpc^^ (Nevalainen & Roos 1998). Since the age of the 
Universe to strongly depends on the Hubble constant, our value also changes 
with Hubble value. Note that the Universe must be older than any other object 
in the Universe. Our determination agrees well with others as can be seen in 



table 1.2. The age determination from the different experiments are summarized 



in table 1.2 



Technique h Assumption Age (Gyr) Objects References 



Stellar ages 
Stellar ages 
High-Z 
High-Z (flat) 
SCP 

SCP (flat) 
4-constraints 
9-constraints 
MCA*) 



None 

None 
0.65 ±0.02 
0.65 ±0.02 

0.63 

0.63 
0.68 ±0.10 
0.68 ± 0.05 
0.70 ±0.15 



13.3 ± 1.1 
12.5 ± 1.5 
14.2 ± 1.7 
15.2 ± 1.7 
14.2 ± 1.0 

14.91}:^ 

13.4 ± 1.6 

13.5 ± 1.3 

1 o n+3.6 
ia.^_2,o 



Halo GC 
Halo GC 
Universe 
Universe 
Universe 
Universe 
Universe 
Universe 
Universe 



m 



m 



100 1 



100 



@ 



Table 1.2: Comparison of the age of Universe estimated using different observa- 
tions. Errors are given at la confidence level. *) Multiband Colour Alalysis of 
bright clusters. 
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Figure 1.1: Geometry of different types of Universe in (O^, ^^a) plane. Flat and 
open models are infinite; closed models are finite. There are also possibilities 
that an infinite Universe can collapse if J^a < 0. A finite Universe expands for 
ever if JIa = with $7^ < 1, and collapses into a big crunch if Clm > 1. A Finite 
Universe can be accelerating/decelerating. 
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Chapter 2 

Problem of Missing Energy 
and Cosmological Constant 



2.1 The Cosmological Constant 

The cosmological constant is not a constant it might be change with time. It is 
related to vacuum energy density, Eq. ( |1.5| ), and it is a potentially important 
contributor to the dynamical history of the Universe. The relation between the 
cosmological constant A and the vacuum energy density were shown in the 
previous Chapter. A strong evidence for a non-zero cosmological constant comes 
from the Supernovae Cosmology Project (SCP), Perlmutter et al. (1999), and 
the High Redshift Supernovae Search Team (HSST), Riess et al. (1998). In 
contrast to standard general relativity, a wide theoretical discussion on a non-zero 
cosmological constant can be found in Carroll & Press (1992), Carroll (2000). We 
show that three independent constraints strongly rule out the standard model of 



flat space with vanishing cosmological constant (Fig. 2^ & Fig, 2^ and for details 
Paper I). However, the present value of the cosmological constant is an empirical 
issue, a precise determination of which would be one of the greatest successes of 
observational cosmology in the near future. If the cosmological constant today 
comprises most of the energy density of the Universe, the age of the Universe is 
much older. 

The most surprising recent advance in cosmology is that 70% of the Universe 
seems to be made of vacuum energy. We have combined most recent observational 
data in the next Chapters to demonstrate this evidence. Our two dimensional 
maximum likelihood gives a matter density of = ^dark + = 0.31 and 
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Figure 2.1: A is in the center of physics. 



a vacuum energy density of f^A = 0.68. The total density of the Universe is 
no = nm + ^A = 0.99 (Paper IV). 

In figure 2^ & Fig.^^ we show our recent (Paper IV) allowed region in 
(r^m, f]A)-plane. The small solid square and the big dashed square are for the flat 
and two-dimensional cases indicated, respectively. 



2.2 The Missing Energy Problem 

The missing energy may be realized by comparing with the critical energy and the 
observed energy. Inflationary cosmology and current cosmic microwave anisotropy 
measurements suggest that the Universe is flat. At the same time observations 
indicate that the sum of ordinary (baryonic) and dark matter in the Universe is 
below the critical density. So, the "missing energy" problem arises. It is possible 
that the missing energy could be from the interacting fields or topological defects 
(Steinhartz 1996, Huey et al. 1999, Bludman & Roos 2001). As the Universe ex- 
pands, the missing energy density varies as {R)^^^^^°'\ where q = for ordinary 
matter, — ^ > a > - 1 for quintessence and a = —1 for the cosmological constant. 
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Figure 2.2: Regions allowed in the (J7oi ^A)-space. The set of solid lines are 
imposed by the observation of the radio galaxy 53W091 assuming a = 0. The 
numbers indicate the value of the Hubble constant in units of km s^^ Mpc~^. 
The limits ilm < 0.22, > 0.6 are indicated. The solid and dashed square 
boxes indicate the recent (Paper IV) allowed region in the parameter space. 
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We can write the luminosity distance relation as 



HodL 



z + -(l-go)^' 



where 



:iim + ( 7, J^^A, 



(2.1) 



(2.2) 



qo is the deceleration parameter, a is the ratio of pressure to energy density of 
missing energy, and its significance is described in Chapter 1. The deviation of 
linear law depends on equation of state a, the ratio of pressure to energy density 



of missing energy as seen in Fig. 2.4 



HodL 




0.4 0.6 

Redshift (z) 

Figure 2.4: The luminosity distance versus redshift in different cosmologies. The 
lower solid line corresponds to Qm = 1, ^^A = and a = — ^ or — |, Einstein-de 
Sitter, the middle solid line corresponds to a = — |, ^Im = 0.31, J7a = 0.68, 
the upper solid line corresponds to a = — |, J7m = 0, Q,a = 1 de Sitter, the 
lower and upper dot lines correspond to a = — i, = 0.31, $7a = 0.68 and 



a = — i, Qm = 0, Oa = 1 de Sitter respectively. The effects of different models 
are observable at high redshift rather than at low redshift. 
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Chapter 3 

Observational Constraints 



While developing science and technology, astrophysical objects become more vis- 
ible to us, leading to remarkable developments in cosmology in the last few years. 
The constraints on cosmological parameters have been published from the anal- 
ysis of different observational data. The constraints on cosmological parameters 
differ from one observation to another. Sometimes the conflict is very strong. We 
partake in resolving this contradiction by paying rigorous attention to statistics, 
and by combining data from recent observations. In a sequence of Papers I-III 
we combined all types of data having a published error on Qm and Oa, using 
simple analysis. This, of course, implied believing in the errors and treating 
them as Gaussian. The techniques and observations usually attributed by the 
astrophysicist for constraining cosmological parameters are summarized in this 
Chapter, and the results of our analysis are discussed. 

3.1 Cosmic Microwave Background Radiation 

The Cosmic Microwave Background (CMB) is one of the important probes of 
the early Universe, and it was dramatically discovered in 1965 by Penzias and 
Wilson. The first result on temperature variations was published by the Cosmic 
Background Explorer (COBE) satellite (Smoot et al. 1992). The analysis was 
done in a two-parameter space: the scalar tilt of density fluctuation of the power 
spectrum Ug, and the CMB quadrupole normalization Q. The power spectrum 
P{k) can be written as P{k) oc A;", where n is the tilt of primordial spectrum. 

The primary aim of the CMB experiments is to determine the power spec- 
trum, Ci, of the CMB as a function of multipole moment I. Given a map of 
the temperature T(h) in each direction h on the sky, the power spectrum can be 
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obtained by expanding in spherical harmonics, 

\aim\ = J ^m(n) T(n), (3.1) 
and then squaring and summing the coefficients, 



If the map covers a patch of the sky that is small enough to be approximated as 
a flat surface, the power spectrum can be written in terms of Fourier coefficients: 

r^= Jdn e-'^^-^Tin), (3.3) 

and then 

Ci ^ {\t(), (3.4) 

where the average is taken over all Fourier coefficients / that have amplitude 
I. Thus, each multipole moment C/ measures the rms temperature fluctuation 
between two points separated by an angle 6 ~ (//200)~^ degrees on the sky. 

For constraining the cosmological parameters, one can fit the data in differ- 
ent parameter spaces. Different groups have analyzed the COBE DMR data. 
Lineweaver (1998) and Tegmark (1999) have analyzed the CMB data up to 
/ ~ 800. In the Cold Dark Matter (CDM) model, Lineweaver (1998) used the 
parameter space as follows: the CDM densities (f^cdm); baryons (J^fe), the tilt of 
scalar fluctuation (n^), the CMB quadrupole normalization for scalar fluctuation 
(Q), the vacuum energy (JIa) and the Hubble parameter (/i). Tegmark (1999) 
has added three more parameters in his analysis, tilt of tensor fluctuations (n^), 
the reionization optical depth (r), and the CMB quadrupole normalization for 
tensor fluctuations (r). Moreover, for the best theoretical model, Tegmark and 
Zaldarriaga (2000) have added one more parameter {i^u) for the neutrino mass 
density parameter in their new analysis. So the parameter space is then ten, and 
the range of all parameters can be seen in Tegmark (1999) and Tegmark and 
Zaldarriaga (2000). However, some of the parameters could be dropped out from 
the analysis because their influence was negligible. The 6-dimensional likelihood 
function was then integrated over the remaining parameters. Overall, the obser- 
vations of CMB are shown by the plot of the multiple moment C/ against the 
multiple I. Recently, de Bernardis et al. (2000) and Balbi et al. (2000) have pub- 
lished the BOOMERANG and MAXIMA-1 data. The angular power spectrum 



from their fitting can be seen in figure 3.1 
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We use Lineweaver (1998) in our least square analysis. The best fit point O 
is located asymmetrically within a wedge at (0.45, 0.35) (Fig. 2 in Lineweaver 

1998) . Denoting the distance from O to an arbitrary point P in the plane by r, 
and the distance from O through P to the la wedge line by ro, our constraint is 
therefore of the form jr^. In the maximum likelihood we fit CMB (Tegmark 

1999) likelihood contours with a fifth order polynomial. The likelihood contours 
in the (Jlmj^A) plane available to us from this compilation correspond to 68% 
(Icr) and 95% (1.96(t) confidence, respectively. 



3.2 Supernovae of Type la Observations 

One of the most energetic phenomena observed in the galaxies is the explosion 
of a star into a supernova. The supernovae are divided into two basic types, 
la and II, and these types are based on the hydrogen spectral lines. Type la 
supernovae are brighter than those of type II. Enormous luminosities of type la 
can be compared with luminosities of a standard candle. The distance-luminosity 
relationship is used to determine the cosmological parameters. 

The key approach to determining the cosmological parameters from the su- 
pernovae observations is in the relation between the apparent magnitude (m), the 
absolute magnitude (M), and the luminosity distance (Riess et al. 1998, 

Perlmutter et al. 1999): 

m{z) = M + 5/o5DL(z;l^„,SlA,i?o) + 25. (3.5) 

One can exclude i^o from Eq. (|3.5D by writing 

m(z) - = 5/oc/Z)L(z;n„,riA) + 25, (3.6) 

where M. = M — 5logHQ. The luminosity distance (Dl) of SN la can be written 

where C is the intrinsic luminosity and is the observed flux of the supernova. 

In the Friedman-Robertson- Walker cosmologies, the luminosity distance at a 
given redshift z is a function of cosmological parameters. Giving constraints on 
these parameters the Hubble parameter Hq, the vacuum energy density Qai and 
the mass density fi^ (Caldwell et al. 1998, Garnavich et al. 1998, Riess et al. 
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1998, Perlmutter et al. 1997, 1999), the luminosity distance is 



Sinn 



DL = cHQ\l + z)\nk 

Qkl^'^ dz {i + zf{i + nmz)- z{2 + z)nA 



Jo 




(3.8) 



where Qk = 1 — Qm — ^Ai sinn is sinh for positive curvature (Tlk ^ 0) and sinn 
is sin for negative curvature {0,^ < 0)- 

The distance moduh /{Hq, ^a) can be obtained from the probabihty 
density function (pdf) of these parameters. The hkehhood for the cosmological 
parameters can be determined from the statistics, where ~ f{Ho, ^m, ^a)- 
The and Ha are independent of the Hubble parameter Hq. Therefore, one 
can convert the three-dimensional pdf to a two-dimensional one as follows 



High-Z Supernovae Search Team (HSST) and Supernovae Cosmology Project 
(SCP) have searched supernovae at different redshift. HSST (Riess et al. 1998) 
has discovered 16 SNe la in the redshift range 0.16 - 0.62 and 34 nearby supernovae 
for the constraining on Hq, Q,m, 0,a, the deceleration parameter (qo), and the 
dynamical age of the Universe (to), whereas the Supernova Cosmology Project 
(Perlmutter et al. 1999) has discovered 42 supernovae in the redshift range 0.18 
-0.83 for the constraining on cosmological parameters. 

The HSST team has used two different methods, a template method and 
Multicolor Light Curve Shape (MLCS) method respectively, for the data fitting. 
The differences of the two methods can be distinguished from the contour plot 
in their fits (their Figs. 6 and 7). A large portion of the contour plot for the 
MLCS method is in the unphysical region $7^ < 0. This can be improved by 
using Feldman and Cousins (1998) statistics. This is a classical statistics, having 
the advantage that is always pulls the fit from the unphysical region to a physical 
region. However, this could only be done when the original data were analyzed. 
Their best value is in the (Q„i,^}/i^) plane (0.20,0.65). The likelihood contours 
in the (0,^, ^a) plane available to us from this compilation correspond to 68.3% 
(Icj), 95.4% {2a), 99.7% {3a) confidence, respectively. We have used this data set 
for both the a-^d the maximum likelihood analyses. In our fit these constraints 
are represented by a term in the x^-sum of the form 



f + OO 



pdf{Q.m,^A)=l pdf{9.m,^A,Ho)dHi 



Jo 



0- 



(3.9) 



[w' + z'-{al-al){wyal)]/al 



z ' 



(3.10) 
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where w and z are the rotated coordinates 



w 

z 



= {flm - 0,mfl) COS 9 + {Qa - ^A,o) sm 9 
= -{flm - ^mfi) sin9 + (JIa - $1a,o) cos 6* . (3-11) 



The rotation angle is ^ = 51°. 1 and the w and z errors are a^j = 1.27, az = 0.18 
for the HSST, whereas for the SCP, the rotation angle is 9 = 54°. 35 and the w 
and z errors are cr^ = 1.3, az = 0.15. The same data set has been analyzed with 
maximum likelihood, and approximating the published confidence contour by a 
fifth order polynomial. 

The sample of the Supernova Cosmology Project is larger than HSST. Due 
to the large sample, the statistical uncertainty is small and the confidence region 
is narrow. It is also possible to estimate some systematic uncertainties in this 
case. In the flat case the value found is Qm = 0.28 it 0.085(stat) ± 0.05(syst). 
The two-dimensional contour plot in the parameter space (0^, ^^a) plane can be 



seen in fig. 3.2. We have used this data in our both analyses. None of the teams 
has considered the effect of galactic dust on the light curve. The effect of inter- 
galactic dust on light curve has been described by Aguirre (1999). However, both 
supernova projects have concluded that the expansion of the Universe is acceler- 
ating. The coming SNAP satellite will be able to measure an enormous amount 
of distant supernovae every year. This can provide us essential information on 
the critical nature of the Universe. 



3.3 Classical Double Radio Galaxies 

Powerful extended classical double radio galaxies are used for constraining global 
cosmological parameters. It is much like the relation between supernovae and 
standard candles, comparing coordinate distance and powerful radio galaxies. 
At first Daly (1994, 1995) applied this method for constraining cosmological 
parameters and later adopted by others Guerra & Daly (1996, 1998), Guerra, 
Daly k Wan (2000). 

It is observed that at a given redshift all radio sources have a similar size. 
This size may be estimated in two different ways: i) by the average size of the full 
population of powerful extended radio galaxies at that redshift; ii) by the product 
of average rate of growth of the source and the total time for the production 
of powerful jets. The two measures depend on the angular size distance to the 
sources, and the ratio of the two measures depend on the cosmological parameters 
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^rru^A) and a model parameter /5. The rate of growth of the radio sources 
strongly depends on the redshift: it increases with redshift, but the average size 
of the full population decreases monotonically with redshift, for a redshift greater 
than 0.5. 

In 1998, Daly, Guerra and Wan studied fourteen samples for constraining cos- 
mological parameters. Guerra, Daly and Wan (2000) have identified 70 classical 
double radio galaxies with a redshift between zero and two. From there they have 
studied 20 of average size for constraining cosmological parameters. However, the 
likelihood contours in the {0,rn,^A) plane available to us from this compilation 
correspond to 68% (la) and 90% (1.64(j) confidence, respectively. These confi- 
dence regions are quite large. This is therefore one of the weak constraints, thus 
affecting very little the total fit. For the analysis we introduced this constraint 
as a mathematical expression like Eq. ( p. 10 ) and Eq. ( |3.11|) , where the rotation 
angle is 6 = 70°. 4 and the w and z errors are cr^ = 0.84, cj^ = 0.33. The same 
data have been introduced in the maximum likelihood analysis as a fifth order 
polynomial. 



3.4 Gravitational Lensing 

According to Einstein's Theory of Relativity light rays are deflected due to the 
gravity of massive bodies. The deflection of light by massive bodies is known 
as Gravitational Lensing. From the gravitational lensing statistics, the matter 
distribution in the Universe can be studied. Distant objects (i.e. QSO's) are 
usually used as a source of light. Galaxies or clusters of Galaxies are used as a 
lensing plane and on the Earth we are the observers. From this geometry one can 
establish a relation between an angular distance and the cosmological parameters 
at a given redshift. 

Different techniques have been applied for constraining the cosmological pa- 
rameters. Using the results from 5 optical quasar surveys, Cheng and Krauss 
(1998) have re-analyzed lensing statistics, and the best fit of Jim gives the max- 
imum likelihood in the range between 0.25-0.55 in the flat Universe. Systematic 
uncertainties is dominant in this analysis, and the uncertainty comes mainly 
from the galaxy luminosity function and dark matter velocity dispersion. From 
the redshift survey of radio and optical data, Falco et al. (1998) have fltted the 
data independently and combined them. The confidence level is so wide that we 
can not even reach the lower limit of the contour. 
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Chiba k, Yoshii (1999) have presented new calculations of gravitational lens 
statistics in view of the recently revised knowledge of the luminosity functions 
of elliptical (E) and lenticular (SO) galaxies and their internal dynamics. They 
applied their revised lens model to a sample of 867 unduplicated QSOs at z > 1 
taken from several optical lens surveys, as well as to 10 radio lenses. In sharp 
contrast to the previous models of lensing statistics that have supported a high- 
density universe with = 1, they concluded that a flat universe with = 
O.S^q'^ casts the best case to explain the results of the observed lens surveys. 
Instead of using the above quoted value, we use the 68% likelihood contour 
in the two-dimensional parameter space of Fig. 8 of Chiba & Yoshii (1999). We 
integrate out the characteristic velocity dispersion (cr*), and we thus obtain the 
one-dimensional 68% confidence range J7a = 0.70 it 0.16. 

A preliminary result from the Cosmic Lens All-Sky Survey (CLASS) has been 
presented by Helbig (2000). Most probably this is the strongest constraint on the 
($7m, JIa) plane from lensing surveys. Although the maximum likelihood is wide, 
this is a strong constraint on Oa. A large portion of the confidence level is in 
the region JIa < which may be considered unphysical. The upper bound of 
the contour may be interesting for Vt\ from this compilation. At this moment 
we get very little benefit from this result. However, we expect that in the near 
future it would be able to provide us with a strong constraint on the cosmological 
parameters. The likelihood contours in the {^m, ^a) plane available to us from 
this survey correspond to 68.3% (Icr), 90% (1.64(7), 95% (1.96c7), 99% (2.58c7) 
confidence, respectively. 

3.5 Cluster Mass Function and Lyo; Forest 

The most common assumption is that the structure we observe today such as 
galaxies, clusters of galaxies and voids results from the growth of primordial 
density fluctuations by gravitational instability. These fluctuations are normally 
assumed to have originated from a Gaussian random process. The overdense 
regions will decelerate faster than the background Universe, because of their 
large gravitational field, resulting in an increase of their contrast relative to the 
background. If this deceleration is large enough, these regions will turn back and 
recollapse on themselves, resulting in the formation of positive density structures 
such as galaxies and clusters. The opposite phenomenon occurs in underdense 
regions. These regions decelerate more slowly than the background Universe, thus 
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getting more underdense, and eventually become the cosmic voids we observe 
today. 

Croft et al. (1998) have presented a theoretical model for recovering the lin- 
ear power spectrum of mass fluctuations from high-redshift Quasi-stellar Object 
(QSO) Lya spectra. They have tested the model by using 16 QSOs. For the 
precise measurement of P{k) Croft ct al. (1999) have combined 3 more sam- 
ples at rcdshift 2.5 with their previous determination. Most recently, toward a 
more precise measurement of matter density, Croft et al. (2000) have used a 
wide number of samples (53 QSOs) of Lya forest spectra in their analysis, dis- 
covered by the Keck telescope at redshifts between 2 and 4. Weinberg et al. 
(1999) estimated by combining the cluster mass function constraint with the 
linear mass power spectrum determined from the Lya data (Croft et al. 1998, 
1999). For = they obtained = 0-46l!lo;io a^nd for a flat universe they 

obtained fJm = 0.34lQ jg. In the flat cosmology Croft et al. (2000) obtained 
O _ f) cn+o.ia 

0.!JU_Q J^Q. 



3.6 Gas fraction in X-ray clusters 

Clusters emit X-rays which indicates that the clusters consist of a large amount 
of hot gas. From the measurement of the gas fraction in X-ray clusters one can 
reach the Universal matter density $7^. It is observed that at a certain radius, 
the clusters are in hydrostatic equilibrium. This radius is known as the virial 
radius. A common deflnition of the virial radius is -R500; outside this radius 
the density drops below 500 in units of the critical density (Navarro, Frenk &; 
White 1995). All matter outside the virial surface are infalling with the cosmic 
mix of components. Then the baryonic mass fraction measured at the virial 
radius must be unbiased. Thus, by measuring the gas fraction near the virial 
radius, one expects to obtain fairly unbiased information on the ratio of 0,m to 
the cosmic baryonic density parameter 0,^. This type of analysis has been done 
before (Evrard 1997) using the best known value of fi^ to derive a value of $7^. 
In Chapter 5, we do the opposite: we use our best estimation of fi^ to determine 
some poorly known parameters. 

For this purpose Evrard (1997) has used a very large sample of clusters: the 
ROSAT compilation of David, Jones & Forman (1995) and the Einstein compi- 
lation of White &; Fabian (1995). He has obtained a realistic value of 

^ (11.8 ± 0.7) . (3.12) 
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This value includes a galaxy mass estimate of 20% of gas mass, and a baryon 
diminution T(500) = 0.85 at i?5oo (a detail definition of T can be found in 
Section 5.4). To get a precise information on the density parameter Qm, we 
need the precise values of Jib and h. The problem lies in the contradiction in 
baryonic density between the Big Bang Nucleosynthesis (BBN) and the Cosmic 
Microwave Background (CMB) observations by the recent balloon experiments 
BOOMERANG (2000) and MAXIMA-1 (2000). However, taking = 0.020 ± 
O.OOl/i^^ (68% CL) from the low primordial deuterium abundance (Buries et 
al. 2001), and h = 0.73 ± 0.07 (68% CL) Gibson k Brook (2001) who have re- 
analyzed Calan-Tololo and the Center for Astrophysics (CfA) SNe la datasets, 
one obtains $7^ = 0.29 it 0.11. 

3.7 Power Spectrum of Galaxy Fluctuations 

Matter in every direction appears to be distributed in high-density peaks sepa- 
rated by voids. The average separation distance is ~ 130/i~^Mpc, which trans- 
lates into a peak in the power spectrum of mass fluctuations. This provides 
a CO- moving scale for measuring cosmological curvatures. Broadhurst & Jaffe 
(1999) used a set of Lyman galaxies at z ~ 3 finding a constraint in the form 
^Irn = 0.20±0.10-|-0.34nA- Roukema Sz Mamon (1999) have carried out a similar 
analysis of quasars, finding Urn = 0.24 ± 0.15 -|- (0.10 it 0.08)J1a. 

The distribution of galaxies is considered to be Gaussian. For the peaks 
measurements, small angle geometry with the normal could provide precise in- 
formation of the co-moving scale. The density parameters Clm, measured 
from these observations lie between the locii of supernovae and cosmic microwave 
background, that is, (341^ — JIa)- 

3.8 Galaxy Peculiar Velocities 

The large-scale peculiar velocities of galaxies correspond via gravity to mass den- 
sity fluctuations about the mean, and depend also on the mean density itself. 
Two catalogs of galaxies have been analyzed for these velocities in order to pro- 
vide information on Clm- the Mark III catalog (Willick et al. 1997) of about 3000 
galaxies within a distance of ~ 70h~^ Mpc, and the SFI catalog (Borgani et al. 
1999) of about 1300 spiral galaxies in a similar volume. Combining the results in 
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these catalogs, Zehavi &: Dekel (1999) quote the constraint 

f^m /i65^ =^ 0.58 ±0.12 , (3.13) 

in the case of flat cosmology, where the error corresponds to a 90% confidence 
level. Taking the index n of the mass-density fluctuation power spectrum to be 
n = 1.0 lb 0.1 (Bond & Jaffe 1998), one obtains the constraint = 0.55 it 0.14, 
where the error corresponds to a 68% confidence level. 

3.9 X-ray Cluster Evolution 

The Clusters of Galaxies are the largest known gravitationally bound and the 
most massive objects in the Universe. Massive clusters are more luminous and 
they emit X-rays. It is observed that the evolution of the number density of rich 
cluster galaxies breaks the degeneracy between the Universal mass density and 
the normalization of the power spectrum as- This evolution is strong with high- 
mass and low-cTg models. The number density of clusters decreases by a factor of 
10^ from z = to z = 0.5. The same clusters show mild evolution in low-mass 
and high-o"8 models, the decrease is by a factor of 10. Bahcall et al. (1997, 1998), 
Eke et al. (1998), Viana k Liddle (1999), Moscardini et al. (2000) have used this 
degeneracy as a powerful tool for the constraint on cosmological parameters. Note 
that this evolution is strongly model dependent. In the next Chapter we have 
analyzed the recent Advanced Satellite for Cosmology and Astrophysics (ASCA) 
data for the study of baryonic dark matter. 

There are still a few more constraints that have been used in the analysis. 
For critical evaluation of the density parameters 0,m and we have done another 
combined study in the Paper (IV). The rest of the constraints have been discussed 
in Paper (IV), and they will be discussed also in the next Chapter. The data 
used in the method are summarized in tables and 

3.10 Data Analysis: Method 

The statistics is one of the strongest mathematical tools that has been suc- 
cessfully used for experimental and observational data analyses. It is used to 
investigate whether the distributions differ from one another, and of course to 
find the best estimation from a set of data. A confidence level is also determined 
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Sources 



r^o = l^m + References 



High-Z SN la 
SCP 

Double Radio Galaxies 
Various CMBR data 
BOOMERANG CMBR 
MAXIMA-1 CMBR 



0.20 ±0.22 0.65 ±0.28 



0.28 
0.09 



+0.19 
-0.13 

+0.29 
-0.09 



correlated 



0.73 



+0.19 
-0.22 



correlated 
correlated 



0.85 ±0.35 



1.01 



+0.26 
-0.25 



1 03+°°^ 
-'-■'-'"^-0.06 

0.90 ±0.08 



@ 



Table 3.1: Two-dimensional data in the analysis. All errors are given at la 
confidence level. 



Sources 



at 



^flat 



References 



Galaxy Power Spectrum 
Gravitational Lensing 
Gravitational Lensing 
Lya Forest 

Temperature- Redshift 
Galaxy Peculiar Velocities 
X-ray Clusters 
Large Scale Structure 
Cluster Evolution 



0.40 ±0.10 
l-f^A 
l-f^A 
"••5^-0.09 

0.27 ±0.10 
0.55 ±0.14 
0.36 ± 0.09 
0.30 ±0.15 
0.45 ±0.20 



1 

0.70 ±0.16 
0.64 ±0.15 

1 
1 
1 
1 

1 ^rri 
1 f^m 



Table 3.2: One-dimensional (flat space) data in the analysis. Errors are given 
at 1(7 confidence level. 



from the combined distribution. The x statistics is calculated as 

X^ = E^^^^. (3-14) 



where Of are the observed and Ej^^ are the expected values. 

Using two free parameters Clrm ^A, we have performed a fit to the con- 



straints summarized in the tables 3.1 and assuming that all reported ob- 
servational errors as well as systematic errors are Gaussian. For this analysis 
we use MINUIT standard program for function minimization and error calcula- 
tion (James & Roos), which is available from the CERN Program Library and 
is documented there as entry D506. We fit the data in two different cases: one- 
dimensional and two-dimensional. The two-dimensional best fit value is found 
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to be = 0.33 ± 0.07, Qa = 0.66 ± 0.12, thus nrn + ^A = 0.99 ± 0.14, where 
the errors are estimated for la (CL). In the exact flatness case we have found 
Qrn = 0.33 ± 0.04, = 0.67 ± 0.04. The interesting thing is that there is no 
remarkable change in the parameters but the errors decrease because the dimen- 
sional space is one. To find the absolute la and 2a confidence regions in the 
Qm, r^A-planc wc add Ax^ = 2.3 and 6.2 up from the minimum, respectively. 
Our fit is extremely good. That can be seen from the degrees of freedom. 
In the two-dimensional fit and for ten constraints, is 4.42, where the degrees 
of freedom is 8, and in the one-dimensional case is 4.43 where the degrees of 
freedom is 9. In the two-dimensional fit is good when is about N - 2, 

where N is the total number of constraints. Therefore the goodness-of-fit here is 
extremely high. 

It is obvious that the greater the discrepancy between observed and adjusted 
values, the greater will be the value of x^. To check this argument we have 
measured the individual x"^ each observational data set. In our ten constraint 
analysis the maximum comes from the CMBR data, but it is still very small: 
slightly less than one. 

We have also fitted the data without the SN la constraints. We found very 
little change. The result is then = 0.34 ± 0.07, Qa = 0.64 ± 0.09, thus 

+ = 0.98 ± 0.11. The reason is that there are strong constraints, for 
instance, gravitational lensing on ^},a and gas fraction in X-ray clusters on the Qfrn 
and they constrain the result similarly. We did not measure any systematic error 
from this fit. If we carefully observe, we do not need to estimate the systematic 
errors. The data we have used here pull in randomly different directions. It 
means that the CMB constraint is laying along the flat line, supernova constraints 
are orthogonal to the CMB constraint. Again the gas fraction in X-ray cluster 
constraints is perpendicular to the gravitational lensing, and the distribution of 
others is quite random. Thus, the total effect might be a mutual cancellation. 
So we can say that the systematic errors arc already included. In such situation 
there is no reason to reconsider the effects of systematic errors. 

There can be more arguments against the systematic errors. Why does one 
neglects them in this analysis? In fact, most of the constraints have come from fits 
of higher dimensions. We do not have any information whether the parameters 
have been properly rescaled or not for the lower dimensionality. Overall, we think 
that our errors have already been generous, and there is no motivation to add 
further systematic errors arbitrarily. 
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Later, we combined a few more constraints in our fit Paper (III), where 
the total number of constraints were sixteen. From the overall analysis, it was 
shown that the flatness of the Universe is robust. 

However, it is observed that a large number of data are not Gaussian, the 
best value is also being unknown in some cases. In such a situation, we loose 
information about the input data. Taking these points into consideration, a 
critical analysis has been done in the Paper IV. The limitations and advantages 
of and maximum likelihood method have been described in Paper IV. Some of 
the constraints that will be used in the maximum likelihood analysis have already 
been discussed in this Chapter. 




Figure 3.1: Top: CMB power spectra from BOOMERANG, MAXIMA-1 and 
COBE-DMR. Bottom: Beam calibration uncertainties are used to adjust the 
peaks (Jaffe et al 2000.) 
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Figure 3.2: The best fit confidence regions in the {0,m, ^a) plane obtained from 
the analysis of the Supernovae Cosmology Project (Perlmutter et al. 1999). 
The 68%, 90%, 95% and 99% confidence region are shown. The diagonal line 
corresponds to flat cosmology. Different cosmologies are also indicated. 
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Chapter 4 

Observational Data Analysis 



Different groups have combined data from different observations, e.g. Lineweaver 
(1999) combined the supernovae data with Cosmic Microwave Background (CMB) 
data, X-ray cluster data, cluster evolution data and double radio sources. Le 
Dour et al. (2000) analysed only Cosmic Microwave Background data, whereas 
Tegmark et al. (2000a) combined them with Infrared Astronomical Satellite 
(IRAS) Large Scale Structure (LSS) data. Tegmark & Zaldarriaga (2000b, 2000c) 
and Hu et al. (2000) combined BOOMERANG and MAXIMA data, Melchiorri et 
al. (2000) combined BOOMERANG and Cosmic Background Explorer (COBE) 
data. Bridle et al. (2000) combined the Cosmic Microwave Background data 
with galaxy peculiar velocities and the supernovae data. The balloon data have 
been combined with LSS and supernovae data by Jaffe et al. (2000) and by Bond 
et al. (2000), and with a different set of LSS data by Novosyadlyj et al. (2000a) 
and Durrer & Novosyadlyj (2000). We have done a critical combined study us- 
ing Maximum Likelihood method. The main interest herein is the following: i) 
taking care of statistics; ii) a combined two-dimensional output contour and the 
best value using individual two-dimensional contours as input; iii) agreement and 
disagreement of data with each other. In this analysis we have used supernovae 
data (Riess et al. 1998, Perlmutter et al. 1999), CMB data (Tegmark 1999, 
de Bernardis et. al. 2000, Balbi et al. 2000, Hanany et al. 2000), LSS data 
(Novosyadlyj et al. 2000b) and Double Radio Galaxies data (Guerra el al. 2000) . 

4.1 Maximum Likelihood Method 

This is another mathematical tool that has been used for the data analysis. From 
the statistical and theoretical point of view, the best method for the data analysis 
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with a high resolution is so far the maximum likelihood method. The hkehhood 
function for a set of data is the joint probabihty density function (pdf) of the 
data, given some parametric model. The values of the parameters maximize 
the sample likelihood. In principle, we can obtain the from the maximum 
likelihood. For a normally distributed variable is the negative logarithm of 
the maximum likelihood. So, we can say that maximum likelihood is optimal, 
and is a special case of maximum likelihood. The advantage of the maximum 
likelihood method is that we can use all kinds of data as an input with more 
precision, details of which may be found in Paper (IV). 

We have represented the confidence contours of the data by fifth order poly- 
nomials of the form 

P(0„J^a) = ^X, (m + n) < 5. (4.1) 

There are 20 terms in the polynomials. Therefore we read off 20 points from 
the la, 1.64o", 2a, 3a contours and the best value, if available. Here we get a 
20 X 20 matrix. We invert the matrix once and multiply with the data points 
that are taken from the contours. The equation is then ready to describe the 
nature of data. Since wc already know the approximate location of the globally 
favored region from all previous studies, we can take care of that our polynomial 
approximation is good over that region. This fit region is defined by 0.15 < 0,jn < 
0.50 and 0.40 < < 0.88, but the sample points are taken also from outside 
this region in order to obtain a well-behaved polynomial inside the region. Far 
away from it, the polynomial approximation of course breaks down completely. 

We have taken special care of each data set. We have checked that the polyno- 
mial is non-negative in the fit region. We have also examined the location of the 
0.33a contour, in order to verify that the region is reasonably centrally located. 
Note that one can use higher order polynomials for better accuracy but it would 
be very difficult to control the polynomial; even in the fifth order polynomial we 
need some constraints to adjust the data points up to the sixth decimal. 

In the previous Chapter, we have described the two well known analyses from 
the supernovae observations and their results. Since we know that the supernovae 
observations constrain essentially only the J7a— parameter combination, let us 
compare their results along the fl^ — 0,^ direction, where they have used the same 
method. Along the flat line the result is then ClA—^m = 0.44ib0.10 for Supernovae 
Cosmology Project (SCP) (Perlmuttcr ct al. 1999), and ^^a - = 0.36 ± 0.10 
for the High Redshift Supernovae Search Team (HSST) (Riess et al. 1998). A 5% 
systematic uncertainty have been included in the SCP result but we do not have 
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any information about the systematic uncertainty for the HSST case. However, 
the two observations then agree within their statistical errors. 

we 



Now we shall see the combined results of these two observations. In Fig. 4.1 



show the confidence contours of the log-likelihood sum of the two observations in 
our polynomial approximation, drawn only in the above mentioned ranges of ^Im 
and r^A. Along the flat line these experiments determine Q\ — = 0.45 it 0.13. 
The SCP result is stronger than that of HSST. Therefore, the best fit is moved 
toward the SCP but the combined contour is naturally wider. 



0/ 




Figure 4.1: Confidence contours of the log-likelihood sum of two SN la observa- 
tions (HSST and SCP). The inner and outer curves corresponds to la and 2a in 
the (Omi f^A)-plane, respectively. The significance of the square is described in 
the text. The diagonal line corresponds to a flat cosmology. 

Constraints on cosmological parameters from earlier Cosmic Microwave Back- 
ground Radiation (CMB) and Double Radio Galaxies has also been described in 
the previous Chapter. By CMB we mean only the observations summarized 
by Lineweaver (1998) and Tegmark (1999), not including BOOMERANG and 
MAXIMA- 1 which we treat separately. CMB puts a strong constraint along the 
direction 0,q = 0,rn + = 1- The combination of supernovae with CMB is also 
well known. It can be interesting to combine Double Radio Galaxies with super- 
novae and CMB. In Fig. [4.2| we show the confidence contours of the log-likelihood 
sum of two supernovae (HSST and SCP), CMB and Double Radio Galaxies in 
our polynomial approximation, drawn only in the ranges of Qm and that we 
sample. The Double Radio Galaxies is so far a weak constraint on and ^}\. 
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Thus the confidence region is mostly dominated by supernovae and CMB. 
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Figure 4.2: The confidence contours of the log-fikehhood sum of two SN la (HSST 
and SCP), CMB and Double Radio Galaxies. The inner and outer curves corre- 
spond to 1(7 and 2(7 in the {^m, fiA)-plane, respectively. The significance of the 
square is described in the text. The diagonal line corresponds to a fiat cosmology. 

4.2 Balloon Experiments 

The recent two balloon experiments BOOMERANG (de Bernardis et al. 2000) 
and MAXIMA-1 (Balbi et al. 2000, Hanany et al. 2000) have given us a most 
exciting information in observational cosmology. The two experiments are held 
in two different parts of the world. BOOMERANG is held in Antarctica and the 
MAXIMA-1 in California. There is a big exposed time difference between two 
observations. BOOMERANG was flown 11 days for the data collection whereas 
MAXIMA-1 was flown only for 2 hours. According to their best estimation, 
BOOMERANG has hinted that the Universe is might even be closed. However, 
MAXIMA-1 tells us that the Universe is open. On one hand we can not escape 
from the statistical fluctuations, and on the other hand there could be some 
important corrections in the theory. Let us now turn to the theory. The widely 
used formula is that the first peak corresponds to a fiat cosmology, and it was 
derived by Kamionkowski et al. (1994) for the case = as 

ii ~ 200^^0"°-^ (4.2) 
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For the case = 0.3 and Oq near 1 Weinberg (2000) has derived the relation 

h ~ 200Jlo~^-^^- (4.3) 

Near = 1 the relation ( [4.2| ) and ( [4.3| ) are of course very similar. Let us now 
turn to their results. BOOMERANG observes (de Bernardis et al. 2000) that 
the position of the first multipole peak occurs at ii = 197 it 6 which corresponds 
to 

= (200/4)^ = 1-01 ± 0.02 . (4.4) 

This value is very weakly dependent on a large number of parameters which 
mostly get determined by the shape of the multipole spectrum above the region 
of the first peak (Lange et al. 2000). Actually BOOMERANG and MAXIMA-1 
fit their multiparameter data by the program CMBfast, so they do not explicitly 



use either (4^) or (O). For the BOOMERANG results, there are no precise 
confidence levels published, only a coarsely pixelized likelihood surface of 95% 
confidence. Other confidence contours such as 68% CL, 90% CL are not reported. 
In consequence we do not have sufficient information to use our polynomial fit 
to the data. Of course the maximum likelihood method can compute these con- 
fidence contours by using the 95% CL information above, but this is not reliable 
at all. Therefore, we did not include the BOOMERANG data as an indepen- 
dent constraint, but it has been included together with the Large Scale Structure 
constraint (Novosyadlyj et al. 2000b). 

Now we turn to MAXIMA-1. From the statistical point of view this is less 
precise than that of BOOMERANG. They find the first acoustic peak at ii ~ 
210. Therefore, we are not surprised about the report of Balbi et al. (2000) 

= 0.90 lb 0.15, where the error corresponds to a 95% confidence level. If we 
convert this to a 68% error, then = 0.90 it 0.08. The likelihood contours in 
the ^m,^A plane from MAXIMA-1 (Balbi et al. 2000) are available to us 68% 
(Icj), 95% (1.96(t) and 99% (2.58(t) confidence levels respectively. 

There is a contradiction in the acoustic peaks between the two observations. 
There is no clear information from where the difference comes. Nevertheless, both 
teams are analysing additional data which may reduce the conflict of the region 
of power spectrum where further peaks are expected. In the meantime, to adjust 
the peaks between these two experiments, a clever technique has been applied 
by Jaffe et al. (2000). Both teams BOOMERANG and MAXIMA-1 have a 
remarkable calibration uncertainty in their method. Jaffe et al. (2000) have used 
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this uncertainty for the adjustment of the peaks. The new acoustic peak is then 
between BOOMERANG and MAXIMA-1. MAXIMA-1 [ and BOOMERANG T 



Fig. 3.1. 



Let us see our confidence region including balloon data. In Fig. |4.3| we show 
the confidence contours of the log-likelihood sum of the SNe la, CMB, Double 
Radio Galaxies and MAXIMA-1. We have not yet included BOOMERANG data 
in this fit. It will be included with LSS in the total fit. We have plotted our 
polynomial approximation only in the chosen window of and 




Figure 4.3: Confidence contours of the log-likelihood sum of two SN la (HSST and 
SCP), CMBR, Double Radio Galaxies and MAXIMA. The inner and outer curves 
correspond to la and 2a in the {^m, f^A)-plane, respectively. The significance 
of the square is described in the text. The diagonal line corresponds to a flat 
cosmology. 



4.3 Large Scale Structure 

Today's great challenge of science is to understand the structure of our Universe. 
One widely used theory on observational cosmology concerns structure formation 
yielding Large Scale Structures (LSS), and it was originated from the theory of 
Big Bang. The form of the power spectrum strongly depends on the cosmological 
parameters. From the observation of the power spectrum of density fluctuations 
one can apply a constraint on Q^n, ^A- In a series of Papers Novosyadlyj et al. 
(1999, 2000a, 2000b, 2000c), Durrer and Novosyadlyj (2000) have analysed up 
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to eight parameters in the fit to obtain the cosmological parameters. They 
have used the constraint on the amphtudes of power spectrum fluctuations from 
different sources on a wide scale. They have combined the power spectrum of 
density fluctuations of Abell-ACO clusters of Retzlaff et al. (1998) with optical 
determinations of the mass function of nearby galaxy clusters by Girardi et al. 
(1998), with evolution of the galaxy cluster X-ray temperature distribution func- 
tion by Viana & Liddle (1999a), Bahcall & Fan (1998), with a study of bulk flows 
of galaxies by Kolatt Sz Dekel (1997), and with Ly-a absorption lines in quasar 
by Gnedin (1998), Croft et al (1998). They apply a constraint on baryon density 
from big bang nucleosynthesis (BBN) by Buries et al. (1999), and the value of 
the Hubble constant h = 0.65 it 0.10. The Hubble value is low but reasonable 
within errors. The position of the first acoustic peak from the BOOMERANG 
observations, £i = 197 it 6, has also been included in their analysis. 

There are seven continuous free parameters in this fit. They have fitted the 
parameters for a fixed neutrino number density (A'^i, = 1, 2, and 3) in a mixed 
dark matter model with a cosmological constant (AMDM). The is minimal 
for one species of massive neutrinos, and the values of mass density and vacuum 
energy density are then = 0.371q^5, = 0-69lo2o- The likelihood contours 
in the Om,f^A plane from LSS (Novosyadlyj et al. 2000c) are available to us 
68.3% (la), 95.4% {2a) and 99.73% {3a) confidence level, respectively. 

In Fig. we show the confidence contours of the log-likelihood sum of 
CMB, MAXIMA-1, LSS and Double Radio Galaxies. We have plotted our poly- 
nomial approximation only in the ranges of and that we sampled. As 
can be clearly seen, the likelihood function contains information mainly on fig, 
but it also gives a rather conspicuous upper limit on the orthogonal combination 
ri/^ — Qm- The narrowness in this fit mainly comes from the CMB and MAXIMA- 
1. 



4.4 Fitting and Errors 

Adding up our polynomial approximations to the confidence contours of all the 
data fitted in this section, results in Fig. 4.5, where we show the location of 
the minimum, the la and 2a contours. From this Figure one can read off the 
following results: 

= 0.3lt°;J2 (4.5) 
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0.2 0.4 0.6 o 0.8 1 

Figure 4.4: The confidence contours of the log-hkeUhood sum of CMBR, 
MAXIMA-1, LSS and Double Radio Galaxies. The inner and outer curves corre- 
spond to la and 2a in the {0,m, r2A)-plane, respectively. The significance of the 
square is described in the text. The diagonal line corresponds to a flat cosmology. 

Oa = 0.68 ±0.12, (4.6) 

or alternatively 

Oo = 0.99 ±0.04 (4.7) 

nA-nm = 0.37t°;20. (4.8) 

Of these results, only the determination of JIq is quite precise and worth 
detailed attention. We can conclude from it that a flat universe with = 1 is 
very likely. 

The error of ±0.04 is mainly statistical, although it does contain systematic 
errors due to mildly discordant experiments, as discussed in the previous section. 
To this we have to add a total systematic error which we evaluate as follows. 

Before discussing the systematic errors let us see the results taking into ac- 
count the new result of lensing statistics Helbig (2000). As we know that there 
is a strong disagreement between lensing and supernovae data, it could be inter- 
esting just to see the combined results. We combine lensing data Helbig (2000) 
with fig. the total results is then Qrn = 0.34lg;o3, Qa = OmtoUt °^ alter- 
natively Qo = 0.97t[!:o4 nA-^m = 0.29to:?8. The best value is excluded by the 
supernovae data at la CL, and it is excluded by lensing data at 97% CL. For this 
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Figure 4.5: All the constraints described in this Chapter combined. The '+' marks 
the best fit: {Qm, ^a) = (0.31,0.68). The inner and outer curves correspond to 
la and 2a in the {^m, J^A)-plane, respectively. The diagonal line corresponds to 
a flat cosmology. 

strong contradiction the nice shape of the contour plot is destroyed; that may be 
seen in the Paper IV fig. 4. 

Perlmutter et al. (1998, 1999) have quoted a total systematic error for 
and JIa along the fiat line of ±0.05. We consider that the same error applies to 

the SN la data of Ricss ct al. (1998), although they could not evaluate it from 
their limited sample of SNe la. Displacing both the SN la contours by ±0.05 
along the flat line, we obtain a very small systematic error in the direction 



There are also two kinds of systematic errors inherent in our method of anal- 
ysis. Firstly, we are reading off the coordinates of the confidence contours of the 
different experiments with some finite precision. We estimate this to be 



Secondly, since we only use 20 points to fit the confidence contours of each 
experiment, there is an arbitrariness in their choice; all we require is that the con- 
fidence contours should be well fitted by whichever polynomial. This polynomial 
arbitrariness results in a systematic error estimated to be 




(4.9) 



A2O0 = 0.027. 



(4.10) 



As^o = 0.01. 



(4.11) 
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The quadratic sum of the errors in Eqs. (4.9), ( 4.10| ), ( [4.11 ) is then 

Atot^o = 0.03. (4.12) 

Thus our final result for 0,q is 

Oo = 0.99 ±0.04 ±0.03, (4.13) 

where the first error is statistical and the second error is systematical. Thus our 
total error is ±0.05. 

Let us now turn to the case of exact flatness, fi^ = 1 — ^A- Along the flat 
line the SN la systematic error is somewhat larger than in the case of Eq. ([4.9|), 



Ail7^'"* = ±0.025. 



(4.14) 



Our result is then 



0^'"* = 0.31 ± 0.04 ± 0.04, 



(4.15) 



where the first error is statistical and the second error is systematical. Thus our 
total error is ±0.055. 

The overall constraints on cosmological parameters in year 2001 are summa- 



rized in table 4.1 



Parameters values 



References 





0.31 ±0.06 




Harun-or-Rashid & Roos 2001 


Qa 


0.68 ± 0.07 




Harun-or-Rashid &: Roos 2001 


Qo 


0.99 ± 0.04 




Harun-or-Rashid & Roos 2001 


to 


13.5 ± 1.3 




Roos & Harun-or-Rashid 2001 


Ho 


73 ±7 




Gibson k Brook 2001 




0.020 ± 0.001 




Buries et al. 2001 (BBN) 




U.UZZ_Q QQ3 




de Bernardis et al. 2001 (CMB) 


a 


-1 < but < - 


1 

3 


Perlmutter et al. 1999 (SNe la) 


a 


-1 < but < - 


1 
2 


Bludman & Roos 2001 (Quintessence) 



Table 4.1: Constraints on cosmological parameters. Errors are given at la con- 
fidence level. The Hubble constant Hq in units of km s~^Mpc~^ and to is in 
Gyr. 
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Chapter 



5 



Search 



for 



Baryonic Dark 
Clusters 



Matter in 



5.1 What is Dark Matter? 



Observations tell that more than 90% matter in the Universe is dark, but it is 
mostly obscure to the physicist and astronomer. In general, any form of matter 
which exists in the Universe in a non-luminous form is called dark matter. Two 
types of dark matter exist, namely baryonic dark matter and non-baryonic dark 
matter. The main baryonic dark matter candidates are Massive Compact Halo 
Objects (MACHOs), white, red and brown dwarfs, planetary objects, neutron 
stars, black holes, dust clouds. The MACHO collaboration (Alcock et al. 1997, 
2000) has surveyed the halos in the Large Magellanic Cloud (LMC) and the 
Small Magellanic Cloud (SMC) at different scales. Searching for microlensing 
events toward LMC Alcock et al. (2000) reported that the most likely MACHOs 
mass is between 0.15Mq and O.9M0. Among known astrophysical objects white 
dwarfs is in this mass range, so they arc good candidates for MACHOs. Most 
recently, a large mimber of dead stars have been detected at the edge of our 
galaxy, the Milky Way (Oppenheimer et al. 2001). Extensions of the Standard 
Model of Particle Physics provide non-baryonic candidates to dark matter. Prom 
the cosmological point of view two main categories of non-baryonic candidates 
have been proposed: Cold Dark Matter (CDM) and Hot Dark Matter (HDM) 
according to whether they were slow or fast moving at the time of the galaxy 
formation. The typical HDM candidates are neutrinos of a few eV, whereas in 
the CDM sector, typical candidates are heavy Dirac or Majorana neutrinos in 
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the GeV-TeV mass range or Weakly Interacting Massive Particles (WIMPs) with 
unknown properties. Supersymmetric theories offer a number of candidates, such 
as neutralinos and axions. 

A good object for the study of baryonic dark matter in the Universe is a cluster 
of galaxies. The size of the clusters are big and the edges are not well defined. At 
large distance from the cluster center their properties approach average universal 
properties. For this analysis Advanced Satellite for Cosmology and Astrophysics 
(ASCA) data will be used from the six X-ray clusters A401, A496, A2029, A2199, 
A2256, A3571. A precise value of can be used in different contexts. Our new 
result of from the Chapters III and IV will be used for the analysis of baryonic 
dark matter in these clusters. Note that this data set has also been analyzed by 
Nevalainen et al. (1999, 2000a, 2000b) for constraining cosmological parameters. 
The measurement of gas fraction in these clusters resembles their analysis. 



5.2 Clusters of Galaxies 

Clusters are composed of baryonic and non-baryonic matter. The baryonic matter 
takes the forms of hot gas emitting X-rays, stellar mass observed in visual light, 
and perhaps invisible baryonic dark matter of unknown composition. Let us 
denote the respective fractions as fgas, fgai, fbdm- It has been observed in a large 
number of clusters that fgas is an increasing function of radius (e.g. White & 
Fabian 1995; Markevitch et al. 1997, 1999; Ettori & Fabian 1998; Nevalainen et 
al. 1999, 2000b), approaching the universal ratio at a large radius, as deduced 
from cluster formation simulations (e.g. White et al. 1993, Frenk et al. 1999, 
Tittley & Couchman 2000). Denoting the universal baryonic density parameter 
as Ofc and the universal total mass density parameter as one can write 

fgasir) + fgalir) + h^r) = T(r)^ , (5.1) 

where T(r) describes the possible local enhancement or diminution of baryon 
matter density in a cluster, compared to the universal baryon density. One can 
make different conclusions from the Eq. f).l[ Let us turn to known and unknown 



quantities in Eq. (5T). The value of is precisely known (Paper IV) and T 
is known with error and there is a strong constraint on it from simulations by 
Frenk et al. (1999). If we trust observational data, fgasir) is known from fig. 5A 



the value of baryonic density Qh is also known but contradictory between BBN 
and CMB (as has been discussed in the previous Chapters). The unknown or 
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poorly known quantity in the above expression is fbdm- Until the contradiction is 
resolved, the baryonic density can be used from BBN and CMB independently. 

It has been known before that the gas fractions in clusters of galaxies are too 
high for Eq. to be satisfied with $7^ = 1 (e.g. White et al. 1993; David et 

al. 1995; White & Fabian 1995; Ettori & Fabian 1998; Mohr et al. 1999). The 
qualitative conclusion given by the lower accepted limit for Q^n is that a large 
fraction of baryonic matter other than gas is ruled out. Using up-to-date data on 
fgas, ^b, ^m, and Hq One can set a quantitative limit to the sum of fgai and 
fbdm using Eq. (O). 



5.3 Gas Fraction 

Usually the total masses and gas fractions have been determined using the hydro- 
static equilibrium condition with an isothermal assumption for the cluster gas. 
Here we have abandoned the isothermality assumption and instead derived the 
gas mass fraction fgasi^) from a sample of six clusters, whose total mass profiles 
have been determined using the gas temperature profiles observed with ASCA. 
For this derivation I am indebted to Jukka Nevalainen. This is an important im- 
provement upon the earlier work, because the accuracy of the total mass within 
a given radius is approximately proportional to the accuracy of the gas tempera- 
ture at that radius. The gas mass profiles have been obtained from the ROSAT 
imaging data (see Nevalainen et al. (2000a) and the references therein for the 
original hydrostatic mass analyses of the individual clusters in the sample). 

We averaged the individual fgasif) profiles without weighing them with their 
errors because the errors are not comparable due to different modeling of the 
temperature data in different cases. As the la error of the average we took 
the standard deviation of their distribution (the square root of the unbiased 
variance of the data set). As a radial parameter we used the overdensity, or 
the mean total mass density within a given radius in units of the critical density, 
d{r) = {p{r)) /p„ where {p{r)) = Mtot{< r)/(|7rr3) and p, = 3F2(l + z)3/87rG {z 
is the redshift of the cluster). Fig. |5]l| shows that fgas increases with decreasing 
overdensity (or increasing radius) reaching a value 

fgas{< rsoo) = (0.200 ± 0.027)/i^o^/' (5.2) 

at a radius where the overdensity is 500. The radial increase is due to the steeper 
decrease of the dark matter density, compared to the gas density in these clusters. 
At the largest radii the fgas values increase rapidly, due to the decline of the gas 
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temperature profiles. In the overdensity range [10^ - 500] the measured fgas 
profile can be approximated by an analytical formula 

fgas{6) = [(V500)-0-0^ - 0.80]h-,^/' , (5.3) 

where the relative error increases from 14% at rsoo to 22% at rioooo- The sys- 
tematic uncertainties inherent in the hydrostatic mass determination method due 
to deviations from the hydrostatic equilibrium or spherical symmetry have been 
evaluated by simulations (e.g. Evrard et al. 1996; Schindler 1996; Roettiger et 
al. 1996). The cluster sample used here has been selected for its lack of any 
signature for such deviations. These uncertainties are negligible compared to the 
above rms variation. 

Our result is consistent with isothermal analyses of cluster samples by Ettori 
& Fabian (1998) and Mohr et al. (1999). For comparison, we also derived the 
average fgas profile for our sample deriving the total masses with hydrostatic 
equilibrium equation, assuming that the gas is isothermal. The resulting profile is 
shown in Fig. 5.1. At small radii the isothermal fgas is higher than the measured 
value, and at large radii the opposite is the case. This is due to the fact that 
at small radii the measured temperatures and consequently the total masses are 
bigger than the isothermal values, but at the large radii the temperature profile 
drops below the average. Even though the isothermal profile is within la errors 
of the measured profile (the isothermal profile gives fgas{< r^m) = (0.170 ± 
0.025)/i5Q ), the data indicates that at radii larger than rsoo the isothermal 
values would probably be significantly smaller than the measured values. 

Due to limitations of ASCA, the measured temperature profiles do not extend 
very far, and consequently the mass determinations are reliable only up to ~ ^500, 
at which radius we evaluate Eq. (p.l|). 



5.4 Local vs. Universal Baryon Fraction 

The cluster formation simulations give information on T(r), or how the cluster 
baryon fractions relate to the universal baryon fraction. Frenk et al. (1999) 
have simulated the formation of an X-ray cluster in a cold dark matter universe 
(assuming Vtm = !> = 0) via hierarchical clustering using 12 different codes. 
Tittley & Couchman (2000), tested the standard cold dark matter cluster forma- 
tion via hierarchical and non-hierarchical clustering. Eke et al. (1998) simulated 
cluster formation in a flat, low density Universe (17^ = 0.3, JIa = 0.7). All these 
simulations give very similar values for T(r) below unity, but approaching it at 
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Figure 5.1: The average fgas profile for the sample with 1 a errors as a function of 
the overdcnsity 6 (solid line). The analytical approximation is shown as a dashed 
line. The profile derived assuming gas isothermality is shown as a dotted line. 

large radii. At the radii corresponding to rsoo in our cluster sample (1.3 - 1.7 h'^^ 
Mpc), the central 68% of the simulated values vary in the range 

T = 0.9 ±0.1. (5.4) 



5.5 Baryonic Dark Matter 

Prom a sample of six clusters with measured temperature profiles we have deter- 
mined an average fgas profile as a function of overdensity. The measured fgas 
values are systematically, but not significantly, below the values determined with 
the usual assumption of gas isothermality at small radii, and above them at large 
radii. 

We combine our result fgasi< r5oo) = 0.200±0.027/i5j'/^ with that of Buries et 
al. (2001) for the baryonic density parameter fife, with the mass density parameter 
^rn value from Paper IV, with the Gibson & Brook (2001) value for the Hubble 
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constant, and the Eke et al. (1998), Frank et al. (1999) and Tittley & Couchman 
(1999) value for T(r). We find the sum of fgai and fbdm is then 

fgai{< rsoo) + fbdm{< r^oo) = -0.010 ± 0.027 . (5.5) 

Since the central value is in the unphysical region, it is more useful to express 
this result as an upper limit. In order to set a meaningful upper limit to the 
sum we must make use of the "unified" approach to classical statistics (Feldman 
&: Cousins 1998), which always produces confidence ranges entirely within the 
physical region. Then we obtain 

fgai{< rsoo) + fbdm{< Tsoo) < 0.027 (84%CL) 
fgai{< rsoo) + fbdm{< Tsoo) < 0.045 (95%CL). (5.6) 

Accepting all the input values, we conclude that the sum of the fractions 
of stellar matter fgai and baryonic dark matter fbdm within rsoo is very small, 
< 0.027, at 84% CL, < 0.045, at 95% CL. Our quantitative conclusion is in 
good agreement with qualitative results in earlier work. In particular, fbdm is 
of the order of a percent, unless the input parameters are wrongly estimated or 
conspire accidentally. But if Qi, were taken to be as small as was estimated by 
Izotov et al. (1999) the situation would be further aggravated. Their estimation 
is ^bhiQo = 0.017 ± 0.003 68% CL), keeping all the values same, then the result 
is fgai{< f5oo) + fbdm{< i^boo) = —0.021 it 0.030. It is also reasonable to conclude 
an upper limit from this result, thus 

fgal{< r-500) + fbdrn{< TSOo) < 0.029 (84%CL) 

fgai{< rsoo) + fbdm{< Tsoo) < 0.048 (95%CL). (5.7) 

Let us use the baryonic density Qb from the recent balloon experiments, de 
Bernardis et al. (2000) and Balbi et al (2000), keeping all the other values same, 
then the result is fgai{< ^^500) + fbdm{< '"soo) = 0.066 it 0.046. We can also 
conclude as a upper limit in the same way, that 

fgai{< r-50o) + fbdm{< ^soo) < 0.049 (84%CL) 
fgaii< rsoo) + fbdmi< Tsoo) < 0.078 (95%CL). (5.8) 

This value is rather high compared with the BBN case. Another discrepancy 
has been found Sadat & Blanchard (2001) between observations and numerical 
simulations concerning fgas- Simulation suggests that the baryonic fractions in 
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X-ray clusters of galaxies are less than the observations found. In this situation 



our values in Eqs. 5.6 and 5.5 will be increased. This discrepancy could be 
resolved from the precise data of X-ray Multi-Mirror Mission (XMM). However, 
it is clear that Qhh'^ needs to be revised. Most recently DASI (Pryke et al. 2001) 
and BOOMERANG (de Bernardis et al. 2001) have reported their flbh'^ values, 
which is a good agreement with BBN, their estimation is Q^/i^ 
68% CL. 



0.0221q;qq3 at 
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Chapter 6 

Black Hole Cosmology 



6.1 Spherical Black Holes 

The theoretical idea of black holes is more than 200 years old, but due to the 
development of science and technology, black holes may become observable and 
interesting nowadays. New information about black holes is given by Chandra, 
who finds a more than 500 solar mass black hole in the M82 galaxy, which is 
located at a distance of 600 light years from the galactic center. This is the first 
confirmed case of such a large black hole outside the center of a galaxy (Chandra 
2000, Kaaret et al. 2000, Matsushita et al. 2000). 

Recent evidence indicates that most of the galaxies have Super-massive Black 
Hole (SMBH) at their centers (Marconi et al. 2000, Peterson & Wandel 2000). 
The formation of a black hole can follow different routes: i) A black hole is formed 
when a collapsed star has more mass than 3 solar masses, ii) A cluster of star- 
like black holes forms and eventually merges into a single black hole, iii) A single 
large gas cloud collapses to form a black hole, iv) A Red giant can explode into 
a supernova and become a black hole. In the case of SMBH, the second case 
is more likely. After the formation of a black hole, it starts to accrete matter 
around it Thorne (1994). It could be interesting to observe the accretion of dark 
matter (DM) in the SMBH. 

We consider that the velocity distribution of dark matter vdm is isotropic. 
In the case of a non-rotating spherical black hole, = Svc, where rj is the 
gravitational radius and Tc is the core radius. The angular momentum of the 
dark matter spiraling at rj is 

jCdm = rriDM vdm n (6.1) 
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where m^M is the mass of dark matter. 

According to classical calculations, the condition that matter can spiral into 
the hole is that the rotational energy of matter should be less than or equal to half 
of the gravitational potential energy (Longair 1994). This condition is fulfilled at 
radius rj 

21 - 2n 

where / is the moment of inertia / = ruDMrf, and the condition that DM falls 
into the black hole is 

^DAi < mDAiVGMn. (6.3) 

The rate of DM falling in the black hole is 

moM = f PDM VDM 47r rf (6.4) 

where / is the fraction of particles falling into the hole. The above equation 
may be written in terms of the universal dark matter density i^oM and the 
Schwarzschild radius. After integrating we get 

Mo - M, = /[144 TT VLdm Pc VDM to]^^ (6.5) 
where Mj is the initial mass , Mq is the present mass, and Iq is the life time 



of the black hole. For /, a simple geometry may be considered, from Fig. B.l 
v± = VDMsinO. Using the angular momentum condition and replacing c = vdm-, 
one may apply a constraint on < 0.42, leading to 

20 

f < — < 0.27 (6.6) 

TT 

Now all parameters are known and one can use Eq. (|6.5| ) to determine the 
accretion of dark matter into a super-massive black hole in a fixed time. For 
instance if a SMBH of mass 10*5 M0 accretes DM for 1 billion years then we get 
from Eq. (U) Mq - Mi k, 10^^ gm. 

According to theory, a black hole emits radiation like a black body. Let us see 
how fast it is evaporating. Collins et al. (1989) have shown that the time needed 
for a black hole to completely evaporate, i.e. to lose all its rest-mass energy M(? 
through such radiation is 

G'^M^ 

t^-f^. (6.7) 
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Figure 6.1: Geometry of a Super-massive Black Hole (SMBH). Particle spiraling 
and the condition that the particle falls into the hole are shown. 



where fi is the Planck constant. The rate of evaporation is 

dM hc^ 



(6.^ 



dt 3G2M2' 

The evaporation depends on the mass of the black hole. So we can conclude that 
the evaporation is more effective in the case of a primordial black hole, and can 
be negligible for the SMBH case. 



6.2 Time of Free Fall 



As we can see from the Fig. 6J dark matter particles are rotating around a black 
hole, thus one can also calculate the rotational time of dark matter particles 
(tDAl) from the above geometry. Dark matter particles rotate around the black 
hole until they reach close to the core of the black hole. Particles gain kinetic 
energy as they approach toward the center, and are finally absorbed by the core. 
Particles can also directly fall at any time into the black hole if they lose the 
angular momentum on their way. 

The gravitational potential energy of a dark matter particle moM at a radius 
is [/ = - GMm,DM and the corresponding force is F = -GMj^mL^ so f = 
_GMtph£^ After integration and substituting the limits we get 

\{rl-rh = -^^. (6.9) 
The gravitational radius ri is larger than the core radius rc, i.e. rf >> rj?. 
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Neglecting the core radius rc we can write 



tDM 



3GM 



(6.10) 



Inserting the values of G, M and rj we can obtain the rotational time of dark 
matter particles around a black hole. It is obvious from the Eq. (6.10) that the 
rotational time of dark matter particle depends on the gravitational radius of the 
black hole. 



6.3 Toroidal Black Holes 

The mechanism of the high-energy jets from Active Galactic Nuclei (AGN) is not 
known. Here and in Paper V we study whether the geometry of a Toroidal Black 
Hole (TBH), in contrast to a Spherical Black Hole (SBH), could properly explain 
the particle dynamics of an AGN. Prom the mass accretion around a TBH, one 
can put a constraint on the life time of AGN activity (Paper V) . A generalization 
of the black hole metrics can be given as follows (Smith & Mann 1997): 

ds^ = -{V + b- — ]dt^ + 



+ R^[d9^ + c sinh^{y/ae)d(f)'^], (6.11) 

where t, r are the time and radial coordinates, and (j) ci-re coordinates on a 
two-surface of constant curvature and F is a potential term. An interpretation 
of the potential would be 

V = (6.12) 

where A is a negative constant in order to provide an anti de Sitter (AdS) space- 
time. This looks like the cosmological constant, which however has a positive 
value (Perlmutter et al. 1997, 1999; Riess et al. 1998; Roos & Harun-or-Rashid 
2000), and which is a few orders of magnitude smaller than needed to act as an 
effective AdS term. The parameters 6, c, a fix the topology of the structure: in 
particular, if 6 = — a = 0(6 = —a following the solution of Einstein field equations 
in empty space) and a — > 0, c = -|-^, then the topology is that of a torus and 
the space-time is asymptotically AdS. We recall that a Schwarzschild metric for 
an asymptotically flat spacetime exhibits 6 = and -|-1 instead of the potential 
term V . Therefore, such a configuration is bounded to the presence of V . Some 
ideas about its origin will be discuss in the next section. 
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6.4 Lifetime of AGN Activity 



From the matter aecretion around a TBH one can obtained the hfe time of AGN 
activity. The AGN activity is related to the toroidal shape of the black hole, so a 
transition to a quiescent state is expected as the black hole reaches SBH status. 
As matter accretes the hole, the event horizon increases, the torus inflates and 
later it looses its starting configuration, turning spherical. On the above line, we 
can infer a more quantitative estimate for the lifetime of the activity phase by 
equating the metric tensor components for a TBH and a SBH. If is an initial 
radial dimension for the torus, e.g. the middle value of the torus thickness with 
respect to the centre of symmetry, and Rfin is the final radius of the SBH, then 
the transition happens if the following condition is matched: 

2MG _ ^ 2MG 

or else: 



V-l^ = ^-l^^ (6-13) 



Rfin = ^_y\j^ ^ (6-14) 



where Vg is the gravitational radius. Besides, the lifetime is related to the accreted 

dt 



matter AM and to the accretion rate by: 



yielding: 



dtilife) oc -^^^Rf,J. (6.16) 

Since ^ is observationally estimated, we could enter the lifetime debate if only 
we knew the potential. Lacking that, we can only study dt{life) as a function of 
some hypothetical functional form for it: 

1. 1/ ~ const; 

2. V(xl/R; 
S. V (X log R; 

where the first case could be regarded as a cosmological vacuum energy density, 
while the other expressions could mimic the background potential of a surround- 
ing axisymmetric galaxy. In particular, the second functional form refers to an 
embedding Newtonian gravitational field; on the contrary, the logarithmic shape 
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is motivated in order to reproduce the flatness of the galactic rotation curve (Bin- 
ney & Tremaine 1987). In this frame, the formation of a TBH can be deduced 
if a protogalaxy develops a sufficiently extended background potential and the 
collapse to a massive black hole is forced to lead to a toroidal configuration. 




Figure 6.2: AGN lifetime for a TBH model with constant accretion rate: {solid 
line) V ~ const; {short-dashed line) V cx l/R; {long-dashed line) V oc log R . 
The normalization has been arbitrarily fixed and r = 

In Fig. |6.2| , the trend of dt{life) is sketched for the above mentioned different 
potentials, assuming a constant accretion rate dM/dt. The \/R potential gives a 
fixed value independent of Rin, while the V ~ const situation gives an increasing 
function, as the log R potential does. If the AGN lifetimes have small dispersion 
and are strictly focused on a single value, e.g. 10^ — 10^ yr (e.g. Cavaliere & 
Vittorini 2000), then the \/R potential seems to give a likely interpretation to it. 

Black hole mass is increasing rather than that of shrinking. In the case of a 
super-massive black hole Hawking radiation is negligible. Toroidal black hole is 
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a nice tool for study AGN phenomena. Prom the mass accretion scenario in a 
Toroidal black hole, we can conclude that the lifetime of AGN activity is finite. 
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Chapter 7 



Summary 

For the sake of minimum loss of experimental information and rigorous attention 
to statistics, both the least square method (x^) and the maximum likelihood 
method have been applied to the analysis of a wide set of astrophysical data. 

The fit which uses 16 independent data constraints, all having published la 
errors, gives the best value Qrn = 0.33±0.07 and il\ = 0.66±0.12, or alternatively 
rio = + = 0.99 ± 0.14. In the exact flatness case Jl^ = 0.33 ± 0.04 and 
Cl\ = 0.67 lb 0.04 (Paper II, III). The goodness of fit of this analysis is extremely 
high, therefore we do not consider separately the systematic uncertainties in this 
fit. 

The maximum likelihood method has been applied to 6 independent data, 
all having plotted contours of several different confidence levels, provides us the 
result = 0.3llgQg and rt\ = 0.68 ib 0.12, or alternatively = + Oa = 
0.99 lb 0.04. Including systematic uncertainties, the exact flatness case yields 
= 0.31 lb OM{stat) ± OM{syst) and = 0.69 ib 0.04(.sy.st) ib OM{stat) 
(Paper IV). Since the results from both analyses are in good statistical agreement 
we have demonstrated that the least-squares method was correct and unbiased, 
and that the neglect of systematic uncertainties was of no importance. 

Thus both and maximum likelihood analyses confirm that the geometry 
of the Universe is fiat, and the Einstein-de Sitter model with Cl\ = 0, = 1 
is strongly ruled out. Any low-density model with J^a = is ruled out. In 
addition, the age of the Universe is measured with a good precision to be to = 
13.5 lb 1.3 {0.68/h) Gyr, whereas to = 12.6 ib 1.2 (0.73//i) Gyr. Ah errors are 
given at la confidence level. 

From an analysis of six clusters of galaxies one concludes that clusters contain 
very little baryonic matter in excess of what is seen in X-rays. Making use of 
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known values for the baryonic density parameter $7^ and previously derived fl^, 
one can establish a limit on the amount of dark baryonic matter in clusters. 
In order to get a meaningful result, the use of a special statistical technique is 
required. 

Black hole masses increase due to the accretion of matter around it. In the 
case of super massive black holes, the evaporation rate is negligible compared to 
the accretion rate. A Toroidal Black Hole (TBH) study in contrast to Spherical 
Black Holes (SBH) shows that the TBH can be used as an important tool in 
explaining AGN phenomena. The acceleration of particles, production of jets, 
the shape of the magnetic field and the lifetime of AGN activity are studied. 

Last few years physicists have been discussed the problem of cosmological 
constant with a great interest but the solution to this problem is still unresolved. 
A new model is needed to resolve this problem. Quintessence model has appeared 
in. However, we have to wait for the success of this model. 

The properties of dark matter in the Universe is mostly unknown, and the 
situation has remained unresolved since last twenty years. The dark matter 
problem can only be partly solved by the most recent discovery of dead stars in 
the Milky Way. A large portion is still unresolved, solution of which would be a 
great success for astrophysics and cosmology. 
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